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1 Introduction 



The interest to ladder type models was raised in a beginning of 90-s (see |T[] for a review) 
in connection with high temperature superconductivity problems in metal oxides. It is 
believed that quasi-one dimensional multi-ladder chains of strongly interacting electrons 
reflect the most important aspects of two dimensional systems and also can reveal some 
properties of the weak coupling between conducting planes. 

Recently there has been considerable interest in the construction of integrable ladder 
type models motivated by the desire to use powerful technique of Algebraic Bethe Ansatz 
(ABA) p|, ^ in the exact investigations of the variety of physical phases of the models. 

In the articles integrable ladder models were constructed by extension of the sym- 
metry algebra, in |Q by defining first the ground state and then formulating a model, 
which has it. The higher conservation laws of an integrable models, which contains next 
to nearest neighbour interactions, were used in construction of ladder models in by de- 
veloping the approach of the article 0.The models with alternating spins were considered 
in P, ^. There are also some other attempts in this area [|10], p!T[ . 

Usually integrable models are homogeneous along the chain, namely, the spectral u 
and model parameters are the same in the product of i?-matrices along the chain. It 
is obvious, that if one considers arbitrary shifts of the spectral parameters by some Zi 
in the monodromy matrix we still have an integrable model. But in order to have a 
local Hamiltonian we need to consider shifts with fixed periodicity n, which causes the 
interaction of spins (or electrons) within an amount of n neighbours, leading to n-ladder 
model. The staggered shift of the spectral parameter was first considered in |]12[ in an 
attempt to construct a relativistic invariant massive Tirring Model in a specific limit of 
the homogeneous XXZ-mode\. 

In recent articles [|l^, Q the authors have proposed a construction of some integrable 
chain models with Z2 grading of the states along the chain, as well as along the time direc- 
tions. Hence, a way of constructing an integrable models with staggered inhomogeneity 
was proposed, motivated by the problems of three dimensional Ising model |]TB| and the 
Hall effect The inhomogeneous anisotropic XX Z Heisenberg chain with staggered 
anisotropy parameter ±A was constructed in |T^, while the isotropic t — J model was 
considered in 



The XX Z and t — J models with only inhomogeneous shift of the spectral parameter 
was previously considered in a chain of articles []18|, [1^, EOl BTI], but these authors have 



not analysed the possibility of the inhomogeneity of the anisotropy parameter A (or other 
model parameters). As a result they have the same intertwining i?-matrix as for ordinary 
homogeneous integrable XX Z and t — J models correspondingly and, therefore, the same 
quantum group structure behind. Contrary to this, in the construction presented in the 
articles and [0 the integrable inhomogeneity appeared not only in shifts of the 
spectral parameter, but also in some structural changes. As a result we got modified 
Yang-Baxter equations (YBE), leading to a quantum group algebraic structure different 
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from the usual slq{2) |T5| . 

As it is shown in the mentioned above articles, due to periodic shift of the spectral 
parameters the Hamiltonians of all this models contain at most next to nearest neighbour 
(NNN) interaction terms and therefore can be represented as integrable two leg ladder 
(zig-zag) models. The interaction between the legs of the ladder is represented in the 
Hamiltonian by a topological type terms written as the product of three neighbour spins 
and the anisotropic antisymmetric tensor as 



As it was shown in |30|, the model of free fermions hopping with inhomogeneous 
parameters, in a case when the rotational invariance is preserved, determines a Peierls 
type mechanism of mass generation, which is based on the breaking of the translational 
invariance for the translations of one lattice spacing. Perhaps the above mentioned 3-spin 
interaction terms are responsible for a gaped phase of the system, as it is discussed in the 
articles [E^] in connection with high temperature superconductivity problems. 

In this article, by use of the solution of the staggered YBE for the general Slq{n) 
quantum group presented in |T^, we pass from spin-1/2 XX Z model to three state 
gl{2\l) supergroup based anisotropic t — J model [|^, |2^, |25|, ^ (or Uq{gl{2\l)) Perk- 



Schultz model) with staggered sign of the anisotropy parameter. Following the technique 



of fermionization of spin models developed in one can see, that anisotropic super- 
symmetric t — J model cam be regarded as fermionized version of the anisotropic spinl 
bosonic Uimin-Lai-Sutherland model [^, 29, l7 . 



In Section 2 we present the solution of the staggered YBE for the Slq{l,2) case and 
calculate the Hamiltonian of the model in a ladder form. 

In Section 3 we apply the technique of the Algebraic Bethe Ansatz (ABA) 0, 0] and 
find the eigenvalues and eigenstates of our Hamiltonian. The Bethe equations [BE) for 
the allowed values of the spectral parameters (momentums) are written down here. 



The Yang— Baxter Equations, their solution and the 
Hamiltonian of the Model 



In this Section we follow the technique developed in the articles |jT3|, We modify 



the anisotropic supersymmetric t — J model (which is the Uq{gl{2\l)) Perk-Schultz model 



2^, ^) in order to construct an integrable model with staggered disposition of the 
sign of the anisotropy parameter A. 

The principal ingredient of integrable models via Bethe Ansatz technique is the R- 
matrix. The i?aj-matrix acts as an intertwining operator on the direct product of the 
auxiliary Va{v) and quantum Vj{u) spaces 

Rajiu,v) : Vaiu) ® V,iv) V,{v) ® Vaiu) . (2) 
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For the anisotropic t — J model the quantum and auxihary spaces are three dimensional 
and corresponding i?-matrix (see for example p4|) can be defined by the following formula 



Raj{z) 



qzRa 



q z 



-1 p-i 



where we have introduced a multiplicative spectral parameter z 



(3) 



q defines the anisotropy of the model and, as it is shown in [25 
i?-matrix Rq is equal to 



= e*". The parameter 
the so called constant 
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(4) 



with X = q — q~^. 

It is very convenient for this model to consider its fermionized version by use of the 
technique developed in the article |2^, 0. Let us introduce two Fermi fields Cj^cr, cf^ with 
spin up and down states a =t, i at each site i of the chain and three basic vectors of the 
corresponding spaces Vi as follows 



-) =|o,i), 



+)=ir,o), 1 0)^1 0,0), 



(5) 



numerated as | 1), | 2), | 3) respectively. Therefore we have Z2 graded quantum Vj{u) and 
auxiliary Va{v) spaces with the following parities 

P(l+))= Ml -)) = !, P(|0)) = 0. 



Consider now Hubbard operators X^, 



0.2 



02) (oi 



defined as in formulas (M) and 



h)'\3\ I, with I a). 



(6) 

I 2) 




|-)(o| 

l+)(0| 
|0)(0| 

(i-^t)4 

c+(l -n;) 

(i-^t)(i- 



(7) 



n 



where n^- = c^TCg- is the particle number operator. The demand that the trace of this 
operator acts on auxiliary and quantum spaces as the identity operator means that the 
double occupancy of the sites by fermions is excluded 



^\l)=XZ\l) = {\-n^n^)\l)=\l), m,l 



1,2,3. 



3 



Now we can write down the formula connecting the fermionic -Raj-operator of the 
anisotropic supersymmetric T — J model with the i^^j-matrix defined by the expression-s 
(I) and (I) 

Ra, = Ra, I Ji) I I (Ji 1= i-iy^'''^'^''\Ra,)Z^X:iXjl. (9) 

It is straightforward to check that this fermionic i?-operator satisfies the YBE in the 



operator form p7 



Let us now consider Z2 graded quantum Vj^p{v) and auxiliary K,(t(w) spaces, p, a = 0, 1. 
In this case we have 4x4 i?- matrices, which act on the direct product of the spaces Va,a{u) 
and Vj^p{v), {a,p = 0, 1), mapping them on the intertwined direct product of Va,a{u) and 
Vj,piv) with the complementary (T = (1 — cr), p = (1 — p) indices 

Raj,ap (U, V) : Va,aiu) ® Vj,p{v) Vj,p{v) (g) VaA^)- (10) 

It is convenient to introduce two transmutation operations Li and L2 with the property 
6^ = 62 = '''d for the quantum and auxiliary spaces correspondingly, and to mark the 
operators Raj,ap as follows 

Raj, 00 = Raj, Raj, 01 = Raj, 

Raj,10 = Rajy Raj,ll = R'aj^ ' (1^) 



The introduction of the Z2 grading in quantum space means, that we have now two 
monodromy matrices Mp,p = 0, 1, which act on the space Vp = H^Li ^j,p mapping it 

on = nf=i Vj,P 

Mp -.Vp^V-p, P = 0,1. (12) 

It is clear now, that the monodromy matrix of the model, which should define the 
partition function, is the product of two monodromy matrices 

M{u) = Mo(m)Mi(m). (13) 

Now, because of the grading in the auxiliary space, we would like to construct the 
monodromy matrices Mq^i as a staggered product of the Raj and R''^j matrices. Let us 
define 

N 

Mo{u) = l[Ra,2j-i{u)K],^{u) 
i=i 

N 

Mi(«) = nKw-i(«)K:25(«). (14) 

where the notation R in general means the different parametrisation of the i?-matrix via 
model (A) and spectral (u) parameters and can be considered as an operation over R with 
property R = R. 
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In order to have a integrable model with commuting transfer matrices 



r(-u) = trQtrQ'M{u) = To{u)ti{u). (15) 
for different spectral parameters 

[t{u), t{v)] = (16) 
it is enough to have the following relations for the ro-('u) = trMcr{u), {a = 0, 1) 

T„{X,u) Ti^a{>^,v) = f„{X,v) fi^a{>^,u), cr = 0, 1 (17) 



It is not hard to see, that in order to ensure the commutativity condition (|TB|) the 
R- and ^-matrices in the expression ( p!4D should fulfill the following two Yang-Baxter 
Equations, which in so called check formalism defined by operator Rij = RijPij, {Pij is 
the permutation operator), has the form 



Ri2{u,v)R23{u)Ri2{v) = R2\{v)Ri2{u)R23{u,v) 
and 



, . Si'-i'-a 



Ru{u,v)R'2T{u)Ri2{v) = R23 {v)R'^2{u)R2z{u,v). (19) 

If the Li and L2 operations are trivial this YBE^s are reducing to ordinary ones. There- 
fore one can take an integrable model with the i?-matrix fulfilling Y BE and try to find 
nontrivial Li and 62 operations. For the Slq{N) case this type of solutions were found 
I3[ from where we will extract the solution for the Slq{2\l) case under consideration. 



m 



The solution of the first set of equations (O) with the nontrivial li operation looks like 



with 
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and the bar operation in (|TH|-|T^) is a multiplicative shift of the multiplicative spectral 
parameter z = by the new model parameter h = e*^ as R{z) = ^(e*") = R{hz~^). 
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This definition for tlie li operation means tliat 



= -iR,)ll iK)ll = -iR,)ll iK)\l = -iR,)\l (22) 

If we define tlie operation L2 on R{z) as 

R''{z) = R''{z-'), (23) 

then easy to see that second set of YBEs will coincide with the first set of YBE-s. 

The exphcit expressions for the matrix elements of R{z) and R''^{z) matrices which 



are fulfilling the YBE-s ( p^p^ ) are the foUowings 



Rll 




= R22 




= q-h-' 


- qz, 


p33 
-^33 


(z) =q h-qz \ 


Ru 




= R22 




= Rlliz) 


= Rlli^] 




-Rll{z) = -Rlliz)=z-z~' 


R21 




= Rll 




= Rlliz) 


= —zX, 


Rlli 


z) = Rlliz) = Rf,{z) = -z-'\ 



(24) 



and 





W)n 


= {R' 




= q 


^ - qz, 






mi 


= {R'' 


mi 


= {R''{ 


z))ll = 


{R^{z)) 




mi 


= {R'' 


mi 


= iR'H 


mi = 


—z\, 




mil 


= {R'' 


mA 


= {R'H 


mi = 


-z-'X 



mi = q ^z-qz \ 

ii = {R^mi = iR'Hm2i = -z+z-' 

(25) 

Therefore we have found an integrable inhomogeneous model and our aim now is to 
write down the expression for the Hamiltonian. But in order to have a local Hamiltonian 
it is necessary to have a point uq, where 

Rijiuo) = const. lij, (26) 

with lij being an identity operator. 

/^From the form of i?-matrix (pO[) one can easily see that the value Uq = (or Zq = 
e™o = 1) is just the needed one. Let us notice that though Rij{u)\u=o = —^hj, we have 
Rij{u)\u=o 7^ iij- Because of this, and as future calculations will show, the Hamiltonian 
of our model contains, together with nearest-neighbour interaction terms, also next-to 
nearest-neighbour interaction terms. 

By definition the Hamiltonian of the model is the logarithmic derivative of the transfer 
matrix at that point 

d\\iT{u) , , , 
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For calculations of logarithmic derivative of the Transfer matrix ( |T5| ) we need to insert 
the linear expansions of the operators R{z), R''^'''^{z), R''^ (z), R'''^{z) around the point u = 



R{u) 
R{u) 
R''{u) 

R'^'Hu) 



-X + uH, H = iqkq + iq~^ R~^ , 
R{h) +uH, H = -iqhRg - iq-^h-^R-\ 
k^{h) + uH'\ H'^ = -iqhR^ - iq-^h~^{R^Y^, 
k^{h-^) + uH'\ H'^ = iqh-^R'^ + iq-^h{R'^Y^, 



- -X-uH 
into the expressions of the Monodromy matrix 

M{u) = k,\{u)Ru{u)kiM ■ ■ ■ k^MRT{n)kUu) ■ ■ ■ 

and extract the linear terms in u from the product. By taking into account that 

k'{h)k^{h-^) = {X^ + Asm^9)i, 
with 9 = —i log h as an additive shift of the spectral parameter, we finally obtain 
- iX{X'^ + 4 sin^ e)Hj_ijj+i = 



(2^ 



(29) 



(30) 



(1 - n,-,)/4(g, h^-'^') + n,-^/o(g, h^-'^') 



(1 -nj)fo{q,h) -n 



h{q,h) + U{q,h-^) 



i2n,-l)feiq,h) 



• [(1 - nj+2,a)fj(mod2)+l{q, h) + (-1)^(1 - nj_i,<j)/j(W2)+l(g, h ^)] 
+ C++1,'^^J> [(^ ~ ^j-l,a)/j(mod2)+l(g, h) + (-l)-'(l - nj+2,a/j(mod2)+l(g, h'^)] ] 



+ 



CM CM 



f2{q.h)+h{q,h-'] 



+ f%{q, h)n2j-in2j 



+ 2S7_i5}5;^.ie-'''=+ (-1)^+1 Si:St^^'^-Uq,h) 

{klni}={j-l,j,j+l} 

+ fsiq, h){nj.inj + UjUj+i) 



(31) 



where fi = 1,2. 

In this expression for the Hamiltonian the spin variables are as in usual t — J model 



s- 
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(32) 



the anisotropic antisymmetric tensors e"'^'^, a,b,c = 1,2,3 



and e 



klm 



k, I, m, = j — 1, j, j + 1 are defined as follows 



2\^a3b 



pafeS 



fliq, h)e 



3af) 



':klm 



— e 



Ikm 



tf\k- 
lf\k-l\ 



l\ 



(33) 



where e is ordinary antisymmetric tensor, and for functions fr{(li h) we got 



/o(g,2: 



2{h 

2{h - h-^){q' - q 
A{h-^ - h){?,qh + qh-^ - g'^/i - 3g"^/i-^) 

2(gl/2/,-l _ ^-1/2^^)2 _ 2(^1/2 _ g-l/2)2^ 

2(g + g-i)(/i-/i-i)2. 



(34) 



It is convenient now to write the Hamiltonian (|3lD in a ladder form represented graph- 
ically as in Figure |l]. 



2j + l 



2j + 3 




s = 



2j 



2j+2 



2j + 4 



Figure 1: Zig-zag ladder chain 



Let us consider the even (2j) and the odd (2j + 1) points of the chain as a sites (j) 
of two different chains labeled by s = and 1 correspondingly. The Fermi fields will be 
marked now as 



0,1 



(35) 



With these notations it is straightforward to obtain from the expression (|3T|) the 
following ladder Hamiltonian 



- + 4 sin^ d)Hj^, = H'. + W. + Hf^ + 



rhS 



aCj,'^,'^ [(1 - %+^,s+l,<?)/4(g, /l"^"^^') + nj+s,s+l,afoiq, h'^^'^^l] } 



- f-11 



'1 — n 



'j + S,S + l 



h{q,h) + U{q,h-') 



'j + S,S + l' 



+ 



(1 - + (-1)^(1 - n,_i+,,,+i„^)/,+i(g, 

(1 - nj+i,,,s)/,+i(g, + (-l)'(l - nj_i+,,,+i,<,/,+i(g, /i^"^)') 



c^M CM f2iq,h) + f2{q,h \ 

7, ^ Jew, ''■j%,i'^i+i,o 



2 

^ahc 



+ 2S'J^^S'j-_|_^^g^j^S'j^_^]^_^e™'' + ( 11 



{fcZm}={(j,s),(j+s,s+l),(j+l,s)} 

+ f^{q,h)nj^snj+s,s+i (36) 

Though this expression is big but one can recognise four clear terms there. The first 
one is the Hamiltonian of anisotropic t — J model for the each of chains (marked as ifj^) . 
The second is anisotropic t — J term for the rungs (i/J„). The third term (//|f ), which was 
discussed in the introduction, is written for the each triangle of the zig-zag, represents the 
interaction between chains and has a topological form. The last term contains spin-spin 
Heisenberg interaction together with hopping of fermions and can be thought as spin-orbit 
interaction. 

At the end of this Section let us write the q ^ 1 limit of this Hamiltonian , which 
corresponds to isotropic case. In this limit it is necessary to rescale the additive spectral 
parameter u and introduce w = log 2/ logg together with new additional model parameter 
6 = log 2/ log g, which should be kept finite. Then we will obtain 



\s+l 



+ 



1 



4 



(2n 



1) 



— nj^snj^s,s+inj+i.s — (— l)*5?T-j_s^j+s,s+i; (37) 

which looks like an ordinary supersymmetric t — J model for each of chains. But the pres- 
ence of the factor 2n — 1 in front of spin-spin interaction means that it is ferromagnetic or 
antiferromagnetic regarding whether the corresponding site in the other chain is occupied 
by fermion or not. 



3 ABA solution of the model 

The ABA solution for the homogeneous Perk-Schultz model was carried out in 
and in a chain of articles pTl for a model with open boundaries. In this section we will 
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use the technique of Algebraic Bethe Ansatz in order to find the eigenstates and 

eigenvalues of the Hamiltonian (^T]) for staggered inhomogeneity. 
At the beginning we need the definition of the L-operator 

^..(^)™ = {m I R^A^) I m') = (-ir("'W"')(i?,,.(^)):^:"'X„", (38) 

which is a 3 X 3 matrix in a horizontal auxiliary space, with the matrix elements acting 
in quantum space. It has the following expression 

(^..)™'(^) = 

(39) 

/ {qz - q-^z-^){l -ni;)ni \ 
+ {z- z-^){l-ni) Xz-^c^ld; -Xz~\l~n^)c'l 

{z - z-^){l - n^)ni 
Xzc-^ci -{-qz + q^^z-^)ni^{l-ni) -Xz-^c^l - n^) 

+ {z ~ z-'){l - n^){l - n^) 

{z-z-^){l -n^)n^ 
-Az(l-n|)c| -Xzc^{l-ni) +{z - z-^)n^\l - n^) + 

\ {q~^z~qz-^){l-n^){l-ni) J 



The matrix elements of the monodromy operator ([T3|) between auxiliary states | k') and 
{k I can be expressed as a product of Lij matrices 

Mo{z)l, = {k I Mo(z) I A;') =L[,l(^)^^Lo2(z)t^---Lo^(z)Jf-\ 

M^{z)l, = {k I M,{z) I k') = L^oT{z)t^Ll^,{-z)tl ■ ■ ■ LUz)lr' (40) 

which in matrix form looks like 




/ As^ll{z) A,,i2(z) 

Ms{zt,= \ As,i2{z) As,22{z) Bs^z) |, s = 0,l, (41) 

where Ag^ah, Bg^a, Cs,a, Dg] (a, 6=1,2) are operators in the the quantum space. 

The graded property of the model is the origin of the following form of the transfer 
matrix (^) 

Ts{z) = -As,^i{z)- As,22{z)+Ds{z), s = 0, 1. (42) 
Now we would like to take the empty fermionic state 

\^)s=\ 0, 0, 0), =1 0)i, I 0)2. ■ ■ ■ I 0)^„ s = 0, 1, (43) 
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as a "test" vacuum and let us demonstrate that the "vacuum" | Q) is indeed eigenstate 
of the transfer matrix (|4^) 

(44) 



The action of Lq^^z) and L''^j^{z) on the fc-th empty state | 0)^ is easy to calculate 



z-z-^ 
z-z 


-h-^z + hz-^ 






-\z-^cl 
1 -A^-ic+ 



T 

q^^z — qz^^ 



(45) 



\hz-^cl 



I 0),, (46) 



h q z — hqz 



-1 



and as we see it has an upper-triangular form. Because the upper triangular matrices 
form a semigroup the action of the monodromy matrix Mq{z)^i on | fi)o also have an 
upper-triangular form directly following from the expression (^) 



( (-h-^z + hz~^)^ 



(47) 



(z — Z Ma 



-h ^z + hz ^) 2 (z — z 



Bo, 



Bo.2 



V 











{h ^q ^z — hqz ^ 
{q^^z — qz^^)^ 



N_ 
I 2 



\Q)o. 



J 



Following the definitions (pT| - p3D of the li and i2 operations one can easily find the 
actions of LqI,'^{z) and -Z^ofc(^) on | 0)^ by use of formulas 



U^'^{z) = L{z~^), U^{z) = U^{z-^] 
and correspondingly the action of Mi{z)^, on | will be 



Mi(^):::' i n), 

( (h-^z-hz-^) 
i-z + z-') 



(45 



(49) 



N_ 
2 



{h-h-hz-^)^{-z + z-^)^ 5i, 



(hq ^z ^ — h ^qz) ^ 
{q~^z~^ — qz)~ 
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Now we see that | f2) is the eigenstate of t(u) with the eigenvalue u^^\z) 



JV 
I 2 . 



-2{-h'^z + hz~^)^{z - z~^)^ + {h'^q ^z - hqz~^)^{q'^z - qz'^ 

-| -\ N -i N 11 1^11 ^ 

-2{h^ z — hz^ )~{—z + z~ )~ + {hq~ z~ — h~ qz)~ {q~ z~ — qz)~ , 



(50) 

which follows from the expressions (p!5|) , (^2|) , (H), (|4^) . Simultaneously we see, that the 
Bs,i{z) and Bs,2{z) (s = 0, 1) operators are creation operators, while the operators Cs{z) 
acts on I Q) as the annihilation operators. This motivates us to consider the states 

\vi,V2,...Vn\F)o = F^--'''Bo,aM)B,^,,{v2)---Bo,aM \^)o, a, = 1,2- (51) 

as n-particle eigenstates of the transfer matrix t{u). is a function of the spectral 

parameters Vj, which should be found later. 

In order to proceed further we need to rewrite the YBEs ([T8|) and ( p!9D in terms of the 
matrix elements of the monodromy matrix 

= (_l)p(^0(pM+P(™0)Mo-(t;)Mit(«)(F)^;,':^,,(^ (52) 

For our purpose, the interesting components of the equations ( ^2)) are the following 
commutation relations between the operators Abc{u) and D{u) with the Bc{v) 

A,Un)BoAv) = ^^^^5l,c'(^)^0,6'a(n) + ^^^i?l,,(«)Ao,ea(t^), (53) 

o[u,v) b{u,v) 

Ao,MB,Av) = -^-^lJ:^BoAv)AiMu) - ^^^Bo4u)A,,ca{v), (54) 

0[U,V) o{u,v) 

D,{u)Bo,a{v) = 777^i?l,a(^)I^0(n)-^5^^i?l,a(u)I^0(^), (55) 

o'[u,v) o'{v,u) 

Do{u)B^M = -777\i?o,a(t^)I^i(«) + ^J^^i?o,a(«)I^i(t;), (56) 
b'[u,v) b'[v,u) 

B,,aiu)Bo,tiv) = ri^'iu,v)B,,a'Bo,b'iu), (57) 
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where we have defined b{u,v), aa{u,v), b'{u,v), a'^{v,u) as follows 

b{u,v) = z{u,v) — z~'^{u,v), 



aa{u,v) 
b'{u, v) 



'\z{u, v), for a = 1 
-Xz~^{u,v), for a = 2 



5 



z(u, v) — z ^(u, v) 



q ^z ^{u, v) — qz{u, v) ' 

Xz(u,v) £ -I 

— for a = 1, 



q ^z{u,v) — qz ^{u,v)'' 

a'a{v,u) = { (59) 
-^^ '("•^) for a = 2 

z{u,v) — qz~'^ (u,v) '' J 

and 

rt:\u,v)= (60) 

/ qz{u,v) — q^^z^^{u,v) \ 

z{u,v) — z^^{u,v) \z^^{u,v) 

\z{u,v) z{u,v) — z^^{u,v) 

y qz{u,v) — q^^z~^{u,v) J 

Here we have used the notation z{u, v) = exp (i[u ~ v]) = j^. 

The first terms in the commutation relations (^3])-(^) are so called "wanted" terms. 
As we will see they are contributing to the eigenvalues of the transfer matrix Ts{u), s = 0, 1. 
The second terms in the (|53|)-(^6|) are so called "unwanted" terms and their contribution 
should be canceled in order to have an eigenstate. 

Using (|53D -(^8|) we can represent the action of the diagonal elements of the monodromy 
matrix as follows 

Di{u) I Vi, ...,Vn I F)i = 

n " 1 N N 

("^^^ri/TT^ — ^[<l'^z~\u) - qz{u)]^[hq-^z-\u) - h-\z{u)]^ \ Vi,...,Vn \ F)q 

n n 

+ Y.^^^t':':.F"'B,,M n (61) 

k=l j=l,j^k 

and 

[Ai^uiu) + Ai,22(m)] I Vi, ...,f„ I F)i = 
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n 

i=i 

n 
i=l 

n n 

k=X i=l,j^k 



(62) 



where 



(63) 



In the equations (|6lD and (^) the first terms appeared due to the "wanted" terms of 
the commutation relations (|55|)-(|5BD while the second terms result from the "unwanted" 
terms. 

The expression ( |5^ for r*^-*^)^^^^^^^ can be rewritten as str[ln{u,Vn)l{n-i)iu,Vn-i)---hiu,vi)], 
with 



( [qz{u,Vk)-q ^{u,Vk)\e\ x -i/- x 2 



(64) 
\ 



v 



\z^u^ ffe)e2 



\qz{u,vi^i - q^z ^(u,Vk)\e\ 
\{z{u, Vk) - z-^{u, Vk))e\ J 



Here the quantum operators act on the k-th space and have (e|^)^ = (5a"^/3* matrix 
representation. 

As it follows from the equations 
eigenstate of Ts{u) 



and (P^D, in order | t>i, ...,Vn \ F) to be an 



Ts{u) I Vi, Vn\ F) = Vs{u] Wl, W„, ) | f 1, ...,Vn\F) 



(65) 



we should demand 

i) the cancellation of the unwanted terms 



bi...b„ f \ ^ \bi...bn 1 pa„...ai _ g 



[(A.)aV.:t - (AOaV.l:]^' 



(66) 



and 



) F should be an eigenvector of the small (nested) transfer matrix r*^^-* (u) 



r^Ki:7M.Vi,...,v^)F"^ = ^A'^ 



u^'>{u;vi,...,Vn)F'''^-''K 



(67) 
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Once this conditions are fulfilled, the expressions for the eigenvalues of Ts{u) become 
iyiiu;vi, ...,Vn) = 

n " 1 JV N 

= (-1)^11777^ -[q-h-\u) - qz{u)]^[hq-^z-\u) - h-^qz{u)]^ 

i=i 

n ^ 

- (-1)^ Wu—^-^'^) + z-\u)]'^[h-^z{u) - hz~\u)]'^v^^\u- Vi, ...,Vr.), (68) 
f \ b{u,Vi) 

with 

uo{u) = ui{u). (69) 

Therefore, as it is obvious from the preceding analysis, in order to know exactly the 
eigenvalue we need to solve an other eigenvalue-eigenstate problem (^) for the transfer 
matrix t^^^ with the reduced amount of degrees of freedom. That is why the all procedure 
is called Nested Bethe Ansatz (NBA). 



In the article ||32[ the computation of the quantities and A^ for the ordinary t — J 
model was demonstrated and the condition of cancellation of the unwanted terms was 
reduced to some equations, which are defining the first set of Bethe equations. It is not 
necessary to repeat this calculations here, but let us mention the main difference, namely 

JV N_ 

(Ai,n + A,22) I ^) = [-ziu) + z-\u)] 2 [h-hiu) - hz-\u)] 2 I n), 

JV JV 

Di\Q)= [-qz{u) + q-^z-\u)] ' [-h-h{u)q + hz-\u)q-^] ' \ Q), (70) 
which leads to the condition 

n. 



T^'\l:.Ci^k;v,,...,Vr.)F''-'^ = -X II [qz{v,,v,)-q-h-\v,,v,)] 
[q~^z~^{vk) - qz{vk)]^[hq-'^z~^{vk) - hr^qz{v]^\^ 



(71) 



The next step of NABA is the diagonalization of the small transfer matrix for 
a chain of length n in order for F to become an eigenvector of it. Following the ABA 
technique we write the YBE 

fkZ'iu,v)M^'Ciu)M^K'iv) = M«:,(.;)M«tM4t'"(^,^), (72) 
where M^^^^ is a monodromy matrix for the "nested" problem, and 
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/ qz{u,v) — q^^z^^{u,v) \ 

Xz{u,v) z{u,v) — z^^{u,v) 

z{u,v) — z~'^{u,v) \z~^{u,v) 

qz{u,v) — q^^z^^{u,v) J 



(73) 



is the r-matrix from (^) in a braid formalism. 
Let us now take 



m™(.) ^ ( ^i;;;*!! i . (T4) 



and the corresponding trace 

r«(M) = (75) 
From the formulas (|72D one can choose following algebraic relations 

A«(n)5«(.) = ^^^"'"\-^"^"^";"^ E(^)(.)AW(n) 



z('U, f ) — 2; -^(m, ti) 



D«(«)5«(.) = ^"f^';^^-^^"(^;"^ i?W(.)DW(.) 

z[u^ v) — z~^{u, V) 



E(^)(u)5«(i;) = 5«(t;)5«(M), (76) 

which are the YBE for the XXZ model with staggered inhomogeneity, defined in the ar- 
ticle [|13| . The first (second) term in the expression (^) called "wanted" (correspondingly 
"unwanted") term. 
Let us take 



|0), 



(1) _ / 1 




n 



|fi)^^^=|0)S^^..|0)« = (9)|0)«. (77) 



k=l 



as a reference state. In order to find the action of nested monodromy matrix M*^^-* (u) on 
the reference state one should act by lk{u) from (|^) on | 0)[.^'*. 
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After simple calculations we obtain 



i=l 

n 
i=l 

Now let us make the following Ansatz for the eigenstates of r^^^ (u) 



(78) 



I v['\ = . . . I ^])«. (79) 

By use of ([7^) the actions of A^^\u) and D^^\u) on the states ([f9|) is given by 



n 



,(1) 



g z[u, VI ' } — qz ^[u, v 
1 z{u,vl^^) - z-\u,v^^'^) fj^ 



m / 



+5:a«s«(.) n 5^^)(-f ) I ^^)^^\ 
fc=i j=ij^k 

A«(«)i.i^...,.«) 



m / _i _!/ (l)x n 

n gz ')-q 'z \u,vl _i _i (i) 
^ nK 177 iiW{u,Vi) - q z [u,Vi)\\ v\\ . . . .v^^J) 

i=i 2;(M,t;) ')-z ^{u,vl ') ^.^1 



+ 5]aW5«(.) n B^^\vf)\nr\ 



^0) 



fc=i 



The first terms in this expressions are the "wanted" terms, while the second terms 
(with A^,^"*, A^^"*) are the "unwanted" terms of the nested problem and should cancel each 



other in the eigenvalue equation (|67D . Hence we can write the eigenvalues 
of r«(M) 



(1) 



n 



qz{u, V. 



[-A«(^)-D«(..)]|t;W^...,^;«) 
' ^ ' ^ ' ^ ^\{[qz{u^v.)-q''z-\u,v.)] 



n 

(1), \{[z{u,Vj)- z-^{u,Vj)] 



q z[u,vl 'j - qz {u,vl 'j -p 
j=i z{u,vl^'^) - z-^{u,vl^^) 



t^iK ■ ■ ■ , v!^^ 



(81) 
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Setting now u = Vk m ( |8TD and comparing it with the condition ([7T| ) we obtain the first 
set of Bethe equations 



n 



qz{vk,vl^^) - q ^{vk,vl^^) [h ^z{vk)-hqz ^{vk)y^q ^z{vk)-qz ^Vk)]'^ 



\ z{vk,Vi^^) - z-^{vk,vl^') 



-h'^z{vk) + hz-^'ivk)] 2 [z{vk) - z-^{vk)] 2 



k = 1, . . . , n. 



^2) 



The equations which are obtained from the cancellation of the unwanted terms 
and M!^ 



(1) 



m _i / (1) (1) n / (1) ^ -1 -u (1) ^ 

k = 1 m 



^3) 



are the second set of Bethe equations. 

Finally we obtain the eigenvalues of ti{u) and tq{u) from (|68[) and (|69| ) as 

Mu; {v}; = (-l)t f[ ^^[q-^-^u) - ^-i^-^H - h-'qz{u)] 

i=i 

+(-1)^ n + - hz-\up ■ 

i=l ^^'^ 

/ -1 -1/ " 

nqz[u,vl ') - q 'z (m,^^- j-rTr / n -i -i/ m 



and 



+n 



-1 / -1/ »^ 



n 



Z\U^ Vj) — Z [U, Vjj 



Mu;{v};W'^}) = Mu;{v};{v^'^}). 



^5) 
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